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Abstract 



We argue that the leading-order subprocess for Higgs-boson production in associ- 
ation with bottom quarks is bb — ► H. This process is an important source of Higgs 
bosons with enhanced Yukawa coupling to bottom quarks. We calculate the correc- 
tions to this cross section at next-to- leading-order in l/ln(mff/77i&) and a s and at 
next-to-next-to-leading order in 1/ In (mff/mj). 
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1 Introduction 



The standard-model Higgs boson has a very weak Yukawa coupling to bottom quarks, 
proportional to m^/v, where v ~ 246 GeV is the vacuum-expectation value of the Higgs 
field. Therefore, the cross section for the production of the standard-model Higgs boson 
in association with bottom quarks is relatively small at the Fermilab Tevatron [jlj and the 
CERN Large Hadron Collider (LHC) [fj], in comparison with other Higgs-boson production 
cross sections. However, if the bottom-quark Yukawa coupling is enhanced, this production 
mechanism can be a significant source of Higgs bosons |2[ . Such an enhancement occurs, for 
example, in a two-Higgs-doublet model for large values of tan/3 = v 2 /vi, where Vi^ are the 
vacuum-expectation values of the two Higgs fields. A value as large as tan/5 w mt/mj, is 
obtained in the simplest version of an SO(10) grand-unified theory. The Higgs boson may 
be detected via its decay to T+T- |, g § or bb § |, 0, § at the LHC and the Tevatron, 
and at the LHC. 

In this paper we calculate the cross section for Higgs-boson production in association with 
bottom quarks at next-to-leading order. We argue that the leading-order subprocess is bb — > 
H, where the initial 6-quark distribution function is calculated |I| [TT1| . We show that the 
subprocess gb — > Hb is a correction to the leading-order subprocess of order 1/ \n(m H /mb), 
and that the subprocess gg — > bbH is a correction of order 1/ In 2 (m H /m b ). We calculate 
both of these corrections, and confirm the calculation performed by two of us ten years ago 
@. Our new contribution to this part of the calculation is a proper understanding of the 
relative order of the different subpro cesses. 

Once we properly identify bb —>■ H as the leading-order subprocess, it is straightforward 
to calculate the order a s correction from the emission of virtual and real gluons. This 
calculation is performed here for the first time. We thus obtain the cross section for Higgs- 
boson production in association with bottom quarks at next-to-leading order in both a s and 
1/ ln(mff/mfe), as well as at next-to- next-to-leading order in 1/ ln(m#/mfc). Our calculation 
is valid for both scalar and pseudoscalar Higgs bosons. 

Our calculation corresponds to the inclusive cross section for Higgs-boson production 
in association with bottom quarks, integrated over the momenta of the b quarks. It is 
appropriate to use our results to normalize the inclusive cross section from a shower Monte 
Carlo program, such as PYTHIA or HERWIG, which uses bb —>■ H as the hard-scattering 
subprocess .Q 

The remainder of this paper is organized as follows. In Sec. ^|we explain how to properly 
count the order of the various contributions to Higgs-boson production in association with 
heavy quarks. In Sec. |3] we perform the calculation of the 1/ ln(m^/mg) and 1/ ln 2 (m#/mQ) 
corrections. In Sec. (| we calculate the a s correction. In Sec. |5] we present numerical results 
at the Tevatron and the LHC. 

2 Counting orders 



1 PYTHIA and HERWIG use backwards evolution of the initial-state b distribution functions to give the 
initial g — > bb splitting. 
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Figure 1: Feynman diagrams for gg — > QQH. 



2.1 1/ ln(mij/mg) correction 

In this section we explain how to count the order of the various contributions to Higgs- 
boson production in association with heavy quarks. This counting is a generalization, to the 



case with two heavy quarks in the initial state, of the counting developed in Ref. |12] for a 
subprocess with one heavy quark in the initial state (qb —> qt) . The underlying concepts for 
the organization of the calculation were developed in Refs. |10|, |TTfl . 

The actual physical subprocess for Higgs-boson production in association with heavy 
quarks is gg — > QQH, shown in Fig. 1. Imagine that the heavy quark is very light compared 
with the Higgs-boson. When the initial gluon splits into a nearly-collinear QQ pair, the 
amplitude is enhanced by the propagator of the internal heavy quark, which is nearly on-shell. 
Integrating over the phase space of the external heavy quark yields a factor of ln(m# /ttlq), so 
the splitting of a gluon into a QQ pair is intrinsically of order a s ln(m#/mQ) (for m# ^> mg). 
Such a splitting occurs twice in each of the first two diagrams of Fig. 1, once in each of the 
next four diagrams,^ and not at all in the final two diagrams. 

Another power of this logarithm appears at every order in perturbation theory, via the 
emission of a collinear gluon from the nearly-on-shell quark propagator. Thus the expan- 
sion parameter is a s ln(miy/mg), and since the logarithm is large, the convergence of the 
expansion is degraded. 

The convergence of the expansion is improved by summing these collinear logarithms to 
all orders in perturbation theory [[H], [11]. This is achieved by introducing a (theoretically- 

2 In the middle four diagrams, one gluon splits into a QQ pair, the other into QQH. Only the former 
gives rise to a factor of ln(mjj/mg). 
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(a) 



Q' Q 

(b) (c) 



Figure 2: Feynman diagrams for (a) the leading-order subprocess QQ — > H; (b) gQ — > HQ 
(there is also an s-channel diagram, not shown); and (c) QQ — ► H, where the heavy-quark 
distribution function Q is given by the perturbative solution to the DGLAP equations. 
Figs, (b) and (c) together constitute the 1 / 'ln(mjf / 'tuq) correction to the leading-order sub- 
process in (a). 



defined) heavy-quark distribution function, Q(x, /i), and using the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) equations to sum the collinear logarithms. The heavy-quark dis- 
tribution function is intrinsically of order a s ln(ra/r/ra<2) since it arises from the splitting of 
a gluon into a nearly-collinear QQ pair ||12|| . 

Once a heavy-quark distribution function is introduced, it changes the way perturbation 
theory is ordered. The leading-order subprocess is QQ — > H, as shown in Fig. 2(a). It 
is intrinsically of order a 2 s In 2 (m^/mg), since each heavy-quark distribution function is of 
order a s In (m/j/mg). (There is also a factor of the heavy-quark Yukawa coupling, which we 
suppress throughout this discussion.) 

Consider next the subprocess gQ —> HQ (and related subprocesses), shown in Fig. 2(b). 
This subprocess gives rise to a factor of ln(m# / tuq) from the region where the gluon splits 
into a nearly-collinear QQ pair. However, this logarithm has been summed into the heavy- 
quark distribution function in Fig. 2(a), so it must be removed. This is achieved by sub- 
tracting the diagram in Fig. 2(c), which corresponds to the subprocess QQ — ► H, but with 
the heavy-quark distribution function given by the perturbative solution to the DGLAP 
equation for a gluon splitting to a QQ pair, 

2tt \m 2 Q JJx y \y J 
where g is the gluon distribution function, and the DGLAP splitting function is given by 

P qg {z) = \[z* + {l-zy\. (2) 

After the cancellation of the logarithm, the sum of the subprocesses in Figs. 2(b) and (c) 
is of order a s times the order of the other heavy-quark distribution function, i.e., of order 
a 2 s ln(mij/mg). This is down by one power of ln(m#/mQ) with respect to the leading-order 
subprocess, Fig. 2(a), so it is a correction of order 1 / \n(mH / tuq) . 

Now consider the subprocess gg — > QQH, shown in Fig. 3(a). This subprocess gives rise to 
a factor of \n(mg/mQ) when either gluon splits into a nearly-collinear QQ pair. Since these 
collinear logarithms have been summed into the heavy-quark distribution functions, they 
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(c) (d) 

Figure 3: Feynman diagrams for (a) gg — > QQH (the complete set of diagrams is shown in 
Fig. 1); (b),(c) Qg — > QH and gQ — > HQ (there are also s-channel diagrams, not shown); 
and (d) QQ — > H, where the heavy-quark distribution function Q is given by the perturbative 
solution to the DGLAP equations. These diagrams together constitute the l/ln 2 (m#/mQ) 
correction to the leading-order subprocess QQ — > H. 



must be subtracted. This is shown in Figs. 3(b)-3(d): each of the two collinear regions must 
be subtracted, but we must also add back the double-collinear region [Fig. 3(d)], which is 
subtracted twice [Figs. 3(b) and (c)]. Once the logarithms have been cancelled, the sum of the 
subprocesses in Fig. 3 is of order a 2 . This is down by two powers of \n(m,H / ttiq) with respect 
to the leading-order subprocess, Fig. 2(a), so it is a correction of order 1/ In 2 (m^/mg). 

Thus we see that Higgs-boson production in association with heavy quarks contains terms 
of relative order unity, 1/ ln(m#/m<2), and 1/ In 2 (m/f/mg), depending on whether the initial 
state contains two, one, or zero heavy quarks, respectively. These are the only powers of 
l/ln(mjj/mQ) that appear, to all orders in perturbation theory (12 . 



2.2 a s correction 

Now consider the correction to the leading-order subprocess, QQ — > H [Fig. 2(a)], from 
virtual- and real-gluon emission, shown in Fig. 4. Since these diagrams contain two heavy 
quarks in the initial state, they are of order In 2 (m^/mg), i.e., down by one power of a s 
from the leading-order subprocess. 

There is a factor of ln(m///mg) associated with the emission of a collinear gluon from a 
heavy quark [Fig. 4(b)], and this is handled in a similar manner to the collinear logarithm 
associated with a gluon splitting to a QQ pair. The collinear logarithm is summed, to all 
orders in perturbation theory, into the heavy-quark distribution function, and the collinear 
region is then explicitly removed by subtracting the subprocess QQ — > H [Fig. 4(c)], with 
the heavy-quark distribution function given by the perturbative solution to the DGLAP 
equation for a gluon radiated from a heavy quark. After the cancellation of the collinear 
logarithms, the sum of Figs. 4(b) and (c) [as well as Fig. 4(a)] is a correction of order a s to 
the leading-order subprocess QQ — > H. 
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(a) (b) 



Q 

(c) 



Figure 4: Feynman diagrams for (a) the virtual-gluon correction to QQ — > H; (b) QQ — > i/g; 
and (c) QQ — > H, where the heavy-quark distribution function Q is given by the perturbative 
solution to the DGLAP equations. These diagrams together constitute the a s correction to 
the leading-order subprocess QQ — > H. 



3 The 1 / \n(mH I ttiq) correction 



The leading-order hadronic cross section for Higgs-boson production in association with 
heavy quarks, the l/ln(ra#/mg) correction, and the 1/ In 2 (tuh/tuq) correction are obtained 
from the equations 

(3) 
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(5) 

where Q is the heavy-quark distribution function, Q is the perturbative heavy-quark distri- 
bution function [Eq. ([[])], g is the gluon distribution function, and they are convolved with 
the various subprocess cross sections in the usual way. (The parton distribution function 
written before the subprocess cross section is from hadron A, the one written after from 
hadron B, and the superscripts on the subprocess cross sections denote the power of a s .) 
This formula implements the discussion in Sec. |2] on the proper way to count the order of the 
contributions to Higgs-boson production in association with heavy quarks.^ One can check 
that the sum of these equations is equivalent to Eq. (5) of Ref. 0, although the proper way 
to count orders was not known at that time. 



3 The next-to- next-to-leading-order formula, Eq. (0), not only subtracts the double-collinear region, as 



described in Sec. 2.1, but also subtracts the single-collinear regions. 
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For the calculation of the l/ln(m#/mg) correction (and also the a s correction), it is 



more convenient to regulate the collinear divergence with dimensional regularization flll 
rather than with a finite heavy-quark mass [10]. The former is accurate up to powers 



of mg/m^, which is small in the region of validity of our calculation, mg <C We 
perform the calculation of the 1/ \n(mH /mo) correction both ways, analytically in the case 
of dimensional regularization and numerically in the case of a finite heavy-quark mass, and 
find agreement. The calculation of the 1/ In 2 (m^/mg) correction is only done numerically, 
using a finite heavy-quark mass. 

We now describe the analytic calculation of the 1/ ln(mH /tuq) correction using dimen- 
sional regularization. The calculation is similar to the QCD correction to the Drell-Yan 
process from initial gluons |]13| ,0 but with the vector current replaced by a scalar current. 

The spin- and color-averaged cross section for the leading-order subprocess QQ — > H is 



a 



QQ->H g y 



7T 771 (//) 2e 
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-6(1 -z) 



(6) 



H 



where z = m? H / s. The calculation is performed in N = 4 — 2e dimensions; /i is the 't Hooft 
mass, which is introduced such that the renormalized Yukawa coupling is dimensionless in 
N dimensions. We use the MS scheme to subtract ultraviolet (and collinear) divergences, so 
m(/i) is the heavy-quark MS mass. 

The spin- and color-averaged cross section for the subprocess gQ — > HQ is 



a s m 2 (fj,) 



'gQ^HQ 



12 
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+ ln 
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(7) 

where z = m 2 H /s. The collinear divergence manifests itself as a pole at e = 0. The 't Hooft 
mass accompanies both the renormalized strong coupling and Yukawa coupling, which are 
defined to be dimensionless in iV dimensions. 

The perturbative heavy-quark distribution function that subtracts the collinear region in 
dimensional regularization is 



dy 



x 



1 



7 + In 4tt 



This is the analogue of Eq. ([!]), in which the collinear divergence is regularized by a finite 
heavy-quark mass. Either distribution function can be used in Eq. (|J), and both yield the 
cross section in the MS scheme. Using dimensional regularization, the first line of Eq. (|j) 
yields the cross section 



a 



i 

gQ^HQ 



a s m (fi) 1 
12 v 2 I 



P qg {z)\n 



m 



H 



:i-z) 



-(l-z)(7z-3) 



(9) 



This is just Eq. (^) with the term proportional to (1/ e— 7+ln Att) removed by renormalization, 
and the limit e — > taken. The remaining terms in the l/ln(m#/mQ) correction [lines 2-4 
of Eq. (H)] yield the same expression. 



'For a pedagogical treatment, see Rcf. jlj 
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The 1/ In 2 (wij/mg) correction, Eq. (H), can also be calculated analytically using dimen- 
sional regularization. However, we find it simpler to perform the calculation numerically, 
using a finite heavy-quark mass. The perturbative heavy-quark distribution function that 
subtracts the collinear region is given by Eq. ([!]), and tuq is kept finite throughout the 
calculation. 



4 The a s correction 

The calculation of the a s correction to QQ —>■ H is also similar to the correction to 
the Drell-Yan process [13, 0. However, there is an additional feature: the ultraviolet 
renormalization of the Yukawa coupling ||15|| . The electroweak coupling is not renormalized 
in the Drell-Yan process due to a Ward identity which cancels the ultraviolet divergence. 

The interference of the one-loop vertex correction in Fig. 4(a) with the tree diagram in 
Fig. 2(a) yields the spin- and color-averaged cross section 
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(10) 

(Cp = 4/3) which includes both the leading-order and next-to- leading-order terms. The 
Yukawa coupling has been renormalized in the MS scheme, which yields the counterterm 
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where 



C = —fM 

5m 
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5m _ a. /I , \ 

— = C F -^3 -- 7 + ln47T . 
m Air Ve / 



(11) 
(12) 



The cross section displays both a collinear (1/e) and an infrared (1/e 2 ) divergence. 

The emission of a real gluon (QQ — > Hg) yields the spin- and color-averaged cross section 



a QQ^Hg 
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(13) 



where z = m 2 H /s, and the "plus" prescription is defined as usual: 



dz [f(z)] + h(z) 



dzf(z)[h(z)-h(l)\ 



(14) 



When combined with the cross section from virtual-gluon emission, Eq. (0), the infrared 
divergences cancel. The collinear divergence is subtracted by constructing the combination 
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(15) 
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which is the analogue, for virtual- and real-gluon emission, of Eq. (H). The perturbative 
heavy-quark distribution function in Eq. (|15D is given by 



dy 
2ir J x y 
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7 + In 47r 



where 



l + z 2 
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(16) 



(17) 



is the DGLAP splitting function for a quark radiating a gluon. The sum of virtual- and 
real-gluon emission, after the subtraction of the collinear divergence, is 
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18) 



which is the sum of Eqs. (|T0|) and fll3|), with the term proportional to (1/e — 7 + ln47r) 
removed by renormalization, and the limit e — > taken. 

Since the derivation of Eq. fllBD involves the removal of both collinear and ultraviolet 
divergences, there are actually two independent scales (/1) present. To make this explicit, 
consider the leading-order running of the heavy-quark mass and the strong coupling |i~5|: 



/ a 8 (/x) V° 
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l + ^o[o; s (//o)/7r]ln(^ 2 /^) ' 

where (3q = (11 — 2iV//3)/4 (see the Appendix for the next-to- leading-order equation). The 
perturbative expansion of Eq. ( |19|) at next-to-leading order is 



m(/i) = m(/i ) 



4+ 



C F — 3 In 
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(21) 



Using Eq. (21) to replace m(/z) with m(nuv) m Eq. (18) yields our final expression for the 
correction from virtual- and real-gluon emission, 
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a. 
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(22) 



where we now distinguish between the renormalization scale (/iuy), associated with the 
running coupling, and the factorization scale (/i), associated with the parton distribution 
functions. 



5 Results and Conclusions 

Our numerical results for Higgs-boson production in association with bottom quarks at 
the Tevatron (y/S = 1.8 and 2 TeV pp) and the LHC (y/S = 14 TeV pp) are presented in 
Table p].f] All cross sections are evaluated with the CTEQ4M parton distribution functions 
p~6|] . The factorization (/i) and renormalization (fiuv) scales are both set equal to m#. The 



running b mass is evolved from an initial value of m^M^) = 4.25 ± 0.15 GeV |l7|]u using 
next-to-leading-order evolution equations (see the Appendix) [T3] .[] The value of the running 
b mass at various values of p>uv are listed in Table [2|. Using m^m^), rather than m^M^) or 
M b , to evaluate the Yukawa coupling decreases the cross section by about 50%. We use the 
standard-model Yukawa coupling, with no enhancement factor, throughout. Although there 



is a lower bound of about 70 GeV on the mass of supersymmetric Higgs bosons |19|], there 
is no strict lower bound on the mass of Higgs bosons in a general two-Higgs-doublet model 
20| , so we include the results for a few smaller masses. 



We show the percentage change in the cross section from the corrections of order 1 / ln(m#/m&) 
1/ In 2 (m#/rri&), and function of the Higgs-boson mass, in Figs. 5 (Tevatron) and 6 

(LHC). We find that the 1/ ln(?7iff/ra&) correction is significant and negative at the Tevatron 
(LHC), ranging from -94% (-86%) for m H = 40 GeV to -35% (-36%) for m H = 1000 GeV. 
The size of this correction is a measure of the validity of our calculation; as it approaches 
approximately —100%, it is no longer justified to regard bb — > H as the leading-order sub- 
process. Our calculation is increasingly justified as one increases the Higgs mass, but the 
l/ln(m#/mb) correction is significant even for m# = 1000 GeV. 

The a s correction is also significant, and happens to be positive, such that it largely 
cancels the 1 / \n{mH / Trib) correction. The a s correction ranges from +56% for m# = 40 
GeV to +78% for m# = 1000 GeV at the Tevatron. The correction increases as the Higgs- 
boson mass approaches the machine energy due to the presence of a large Sudakov logarithm 



21|.P| Such an effect is not present at the much higher-energy LHC, where the correction 
ranges from +52% for m H = 40 GeV to +42% for m H = 1000 GeV. 

The 1/ ln 2 (mij/mb) correction is modest, ranging from +20% (+17%) for m# = 40 GeV 
to +3% (+4%) for m H = 1000 GeV at the Tevatron (LHC). This supports our counting 

5 The contribution from qq — > bbH is negligible at both machines. 

6 Mf, is the pole mass, which equals 4.64 GeV at one loop. This is the appropriate value to use in m(,(-M{,), 
since the MS mass given in Ref. |l7|] is obtained from the pole mass at one loop. 

7 When the b mass appears in the kinematics or the perturbative heavy-quark distribution function 
[Eq. (0)1, the CTEQ mass (5.0 GeV) should be used. 

8 We do not attempt to sum the Sudakov logarithm. 
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200 400 600 800 1000 
m H (GeV) 

Figure 5: Percentage change in the cross section for Higgs-boson production in association 
with bottom quarks from the corrections of order l/ln(m H /m b ), l/\n 2 (m H /m b ), and a s , as 
a function of the Higgs-boson mass, at the Tevatron. The next-to-leading-order (NLO) cross 
section is the sum of the leading-order cross section and the corrections of order 1/ ln(m#/mft) 
and a s . 

of the order of the various corrections. Since we have not calculated the other next-to- 
next-to-leading-order corrections, of order a 2 and a s x 1/ ln(m#/m{,), we do not include the 
1/ In 2 (m^/mb) correction in our final result, given in Table [I]. 

The largest uncertainty in the cross section comes from varying the factorization scale, 
/i. We show in Figs. 7 (Tevatron) and 8 (LHC) the percentage change in the cross section 
from its central value (fi = due to varying [i between and 2m#, while keeping 

= itlh fixed. The scale variation is generally larger at next-to-leading order than it 
is at leading order ,^ which indicates that the leading-order scale variation is not a reliable 
estimate of the uncertainty. The next-to-leading-order uncertainty ranges from about ±30% 
for m H = 100 GeV to about ±8.5% (±13%) for m H = 1000 GeV at the Tevatron (LHC). 

There is a much smaller uncertainty in the cross section from varying the renormalization 
scale, [luv, between m#/2 and 2m h- The next-to-leading-order uncertainty ranges from 
about ±2.5% for m H = 100 GeV to about ±6% (±4%) for m H = 1000 GeV at the Tevatron 
(LHC). This is significantly less than the leading-order uncertainty of about ±11% at both 
machines. There is also an uncertainty in the cross section of about ±8% from the uncertainty 

9 In these figures, the leading-order cross section is calculated with leading-order parton distribution 
functions Jl6[ . 
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Figure 6: Same as Fig. 5, but at the LHC. 



in the 6-quark mass. 

Since the 6-quark distribution function arises from the gluon distribution function, an 
additional source of uncertainty is from the gluon-gluon luminosity, which depends oiitk 
m 2 H /S [T22|. We use the uncertainty advocated in Ref. J22|: ±10% (y/r < 0.1); ±20% 
(0.1 < y/f < 0.2); ±30% (0.2 < < 0.3); ±60% (0.3 < y/T < 0.4). We combine all 
four sources of uncertainty discussed above in quadrature, and report the uncertainty in the 
next-to-leading-order cross sections in Table |]. 

Our calculation might also be applied to Higgs-boson production in association with top 
quarks P, |2"1 |. However, it is only valid for m# 3> m t , where ti — > H can be regarded as 
leading order and gt — > Ht can be regarded as a small correction of order 1/ ln(m#/m 4 ). For 
rriH ~ rn t , one must regard gg — ► ttH as the leading-order subprocess (along with qq — > ttH). 
The a s correction to these subprocesses has not yet been calculated.^ 

In summary, we have calculated the cross section for Higgs-boson production in associa- 
tion with bottom quarks at next-to-leading-order in l/ln(m#/mb) and a s , and at next-to- 
next-to-leading order in 1/ ln(ra#/mi,). The most important effect of the next-to- leading- 
order corrections is taken into account by evaluating the bottom-quark Yukawa coupling 
using the MS mass evaluated at the Higgs mass. The l/ln(m#/m&) and a s corrections are 
both large, but have the opposite sign, such that the total next-to-leading-order correction 
is relatively modest. 

10 However, it has been calculated in the opposite limit to the one we are considering, namely ran <C m t 
(with m H ,m t < v 7 !) |p5fl . 
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Figure 7: Uncertainty in the cross section for Higgs-boson production in association with 
bottom quarks at the Tevatron, obtained by varying the factorization scale /i about its 
central value, /i = ran, from m H /2 to 2m H . 




Figure 8: Same as Fig. 7, but at the LHC. 
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Table 1: Leading-order and next-to- leading-order cross sections (pb) for Higgs-boson pro- 
duction in association with bottom quarks at the Tevatron (\^S = 1.8 and 2 TeV pp) and 
the LHC (v^S = 14 TeV pp). The next-to-leading-order (NLO) cross section is the sum 
of the leading-order cross section and the corrections of order l/k^m^/mb) and a s . All 
calculations are performed in the MS scheme using CTEQ4M parton distribution functions 
with jj, = [ijjv = m H- The method to obtain the uncertainty in the next-to- leading-order 
cross section is described in the text. 



VS = 1.8 TeV pp ' " VS = 2TeVpp = 14 TeV pp 



m H (GeV) LP NLO (pb) LQ NLO (pb) LP NLO (pb) 



40 


12.30 


7.45 ±4.17 


xlO" 1 


15.10 


9.36 ±5.19 


xlO" 


-i 


3.16 


2.10 ± 1.03 


xlO 1 


60 


2.71 


1.99 ±0.85 


xlO" 


-i 


3.41 


2.50 ± 1.07 


xlO' 


l 


9.94 


7.52 ±2.98 


xl0° 


80 


8.23 


6.52 ±2.41 


xlO" 


-2 


10.70 


8.52 ±3.13 


xlO" 


-2 


4.12 


3.33 ± 1.16 


xl0° 


100 


3.05 


2.56 ±0.86 


xlO" 


-2 


4.05 


3.38 ± 1.12 


xlO" 


-2 


2.01 


1.69 ±0.54 


xl0° 


120 


1.28 


1.11 ±0.35 


xlO" 


-2 


1.74 


1.51 ±0.47 


xlO" 


-2 


11.00 


9.55 ±2.92 


xlO" 1 


140 


5.90 


5.28 ± 1.57 


xlO" 


-3 


8.22 


7.33 ±2.17 


xlO" 


-3 


6.47 


5.76 ± 1.68 


xlO" 1 


160 


2.91 


2.67 ±0.75 


xlO" 


-3 


4.15 


3.79 ± 1.07 


xlO" 


-3 


4.02 


3.61 ± 1.01 


xlO" 1 


180 


1.51 


1.42 ±0.39 


xlO" 


-3 


2.21 


2.06 ±0.56 


xlO" 


-3 


2.63 


2.39 ±0.65 


xlO" 1 


200 


8.21 


7.83 ± 2.49 


xlO" 


-4 


1.23 


1.17 ± 0.31 


xlO" 


-3 


1.78 


1.65 ±0.43 


xlO" 1 


250 


2.04 


2.03 ±0.61 


xlO" 


-4 


3.24 


3.19 ±0.97 


xlO" 


-4 


7.66 


7.23 ± 1.79 


xl0~ 2 


300 


5.84 


6.01 ± 1.75 


xlO" 


-5 


0.98 


1.00 ±0.29 


xlO" 


-4 


3.74 


3.58 ±0.85 


xl0~ 2 


400 


6.10 


6.67 ±2.36 


xlO" 


-6 


1.17 


1.25 ±0.35 


xlO" 


-5 


1.15 


1.14 ±0.26 


xlO" 2 


500 


7.70 


8.86 ±3.05 


xlO" 


-7 


1.70 


1.93 ±0.67 


xlO" 


-6 


4.34 


4.34 ±0.92 


xl0~ 3 


600 


1.07 


1.29 ±0.80 


xlO" 


-7 


2.75 


3.25 ± 1.10 


xlO" 


-7 


1.90 


1.93 ±0.40 


xl0~ 3 


700 


1.54 


1.97 ± 1.21 


xlO" 


-8 


4.72 


5.83 ±3.60 


xlO" 


-8 


9.21 


9.44 ± 1.89 


xlO" 4 


800 


2.28 


3.08 ± 1.90 


xlO" 


-9 


0.83 


1.08 ±0.66 


xlO" 


-8 


4.79 


4.97 ±0.98 


xl0~ 4 


900 


3.42 


4.87 ±3.00 


xlO" 


-10 


1.49 


2.01 ± 1.24 


xlO" 


-9 


2.64 


2.77 ±0.54 


xl0~ 4 


1000 


5.18 


7.82 ±4.81 


xlO" 


-11 


2.70 


3.85 ±2.37 


xlO" 


-10 


1.52 


1.60 ±0.31 


xlO" 4 
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Table 2: The two-loop running b mass at various values of the renormalization scale /jLuv, 
given the initial value m b (M b ) = 4.25 GeV. 



fiuv (GeV) 


m b (iiuv) ( 


40 


3.23 


60 


3.11 


80 


3.04 


100 


2.98 


120 


2.94 


140 


2.90 


160 


2.87 


180 


2.85 


200 


2.82 


250 


2.77 


300 


2.73 


400 


2.68 


500 


2.63 


600 


2.60 


700 


2.57 


800 


2.55 


900 


2.53 


1000 


2.51 



eV) 
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Appendix 

The next-to-leading-order running of the heavy-quark MS mass is described by the equation 



where 



ajjj?)\ ^ [1 + Aia a (//)/7r] 

f 7 7T— A J w 1 , 23 

^ = "^ + T. < 24 > 
Po Po 



00 













i( 


11 - 


¥•) 


, 7o 


= 1 , 


1 








1 




(l02 




, 7i 




16 




~~ 16 



where Nf is the number of quark flavors of mass less than \i (// > /io). For fi < m t , Nf = 5; 
for fi > m t , Nf = 6, with the running mass continuous at fi = m t . 

The relation between the pole mass and the MS mass at next-to-leading order is 

M Q = m Q {M Q )(l + C F ( ^ . (25) 

The next-to-leading-order evolution of the strong coupling is supplied as a subroutine by 
CTEQ ||16[| . An analytic expression can be found in Ref. [17]. 
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